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Abstract. We investigate how the dynamic polarization of the carbon atoms valence 
electrons affects the spatial distributions of protons channeled in the (11, 9) single- 
wall carbon nanotubes placed in vacuum and embedded in various dielectric media. 
The initial proton speed is varied between 3 and 8 a.u., corresponding to the energies 
between 0.223 and 1.59 MeV, respectively, while the nanotube length is varied between 
0.1 and 0.8 /im. The spatial distributions of channeled protons are generated using 
a computer simulation method, which includes the numerical solving of the proton 
equations of motion in the transverse plane. We show that the dynamic polarization 
effect can strongly affect the rainbow maxima in the spatial distributions, so as to 
increase the proton flux at the distances from the nanotube wall of the order of a 
few tenths of a nanometer at the expense of the flux at the nanotube center. While 
our findings are connected to the possible applications of nanosizcd ion beams created 
with the nanotubes embedded in various dielectric media for biomedical research and 
in materials modification, they also open the prospects of applying ion channeling for 
detecting and locating the atoms and molecules intercalated inside the nanotubes. 
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1. Introduction 

Theoretical modeling of ion channeling through carbon nanotubes has reached a mature 
level [H [2l E] whereas experimental realization of this process is still at the preliminary 
stage [H |5]. Nevertheless, theoreticians continue to explore possible applications of 
ion channeling through the nanotubes, beginning with those based on the angular 
distributions of channeled ions, e.g., for the purpose of deflecting ion beams in 
accelerators [HI [7] or determining some structural details of the short nanotubes using 
the rainbow effect |Hl El UHl [H] • Besides, theoretical studies of the spatial distributions of 
channeled ions have demonstrated the possibility of creating nanosized ion beams, which 
could find interesting applications in biomedical research and for materials modification 
[HI [121 [13 [in [IS]. However, an additional application of ion channeling through the 
nanotubes based on the spatial distributions remains to be investigated in the context of 
extending the classical use of ion channeling for materials analysis [H] to the nanotube 
based materials. 

Namely, one of the most powerful applications of ion channeling through single 
crystals was developed to detect and locate their defects. It is based on measuring ion 
dechanneling on such defects in combination with their Rutherford backscattering |16j . 
It was crucial for this application of ion channeling to develop the adequate theoretical 
models for calculating the spatial distributions of channeled ions, and thus enable the 
precise control and interpretation of the measurements, providing the information on 
such defects [I6]. On the other hand, the nanotubes can exhibit numerous defects, 
appearing during their synthesis or due to external agents such as particle irradiation 
[IT] . Moreover, one of the main directions within the nanotube research is related to the 
phenomena of intercalation of the atomic and molecular species in them [IB]. Therefore, 
it is reasonable to explore whether ion channeling can be also used to analyze defects in 
the nanotubes, in particular, to detect and locate the atoms and molecules intercalated 
in them. In a first attempt to address this problem, we present here a detailed study 
of the spatial distributions of protons channeled through the nanotubes, with special 
emphasis on the rainbow effect in the transverse plane with its possible applications. 

Having in mind that the ion channeling regime is most suitable for probing the 
atoms intercalated in the nanotubes involves light ions in the MeV energy range 
|16j . one must emphasize that, for such projectiles, the dynamic polarization of the 
nanotube atoms valence electrons can give rise to the strong image force, which can 
pull the channeled ions towards the nanotube wall [19]. Of course, one expects that the 
intercalated atoms would be placed right there, close to the nanotube wall, in some state 
of adsorption on the nanotube atoms surface. While the image force has not been found 
to play any significant role in channels of single crystals [16], its role has been identified 
clearly in ion-surface scattering [20], and in ion transmission through capillaries in solids 
[2T| [22] . In addition, although the image force plays the minor roles in ion channeling 
through the nanotubes in the GeV and keV energy ranges [H [21 [6], [3, [23], it has been 
shown to influence strongly ion trajectories in the MeV energy range. For example, the 
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image force gives rise to the rainbow effect in angular distributions of protons clianneled 
tfirougli the short chiral single- wall and double wall nanotubes [2^ 125]. 

It has been established that, when the nanotubes are grown in a dielectric medium, 
one can achieve a very high degree of their ordering and straightening. This makes 
such composite structures very suitable candidates to be used in experiments of ion 
channeling through the nanotubes. Therefore, it is important to study the effects of 
the surrounding dielectric medium on ion channeling through the nanotubes in the 
MeV energy range. The materials of interest are AI2O3 [1], Si02 [SHI ETj as well as 
Ni [28] and Pt [5]. In this context, it has been shown recently that the image force 
and, consequently, the rainbow effect in the angular distributions of protons channeled 
through the nanotubes can be strongly modified by the polarization of the surrounding 
cylindrical dielectric boundary [29l |30] . The separation between the nanotube and the 
Si02 and AI2O3 boundaries is approximated by the nanotube atom van der Waals radius 
[3T| [32] . The separation between the nanotube and the Ni boundary is approximated 
by the separation between the a graphene and a (111) Ni surface, which was calculated 
using a minimization structural calculation based on the density functional theory |33j . 
In both cases the obtained separation was 3.21 a.u. 

We shall investigate here how the spatial distributions of protons channeled in 
the (11, 9) single- wall nanotubes placed in vacuum or embedded in various dielectric 
media are affected by the dynamic polarization effect. While such a study complements 
our previous studies of the influence of the image force on the rainbow effect in the 
angular distributions of protons channeled in the nanotubes [211 ESI SO] , we expect that 
the analogous influence on the spatial distribution of channeled protons may lead to 
a significant change of the proton flux within the nanotube, enabling one to probe 
the atoms and molecules adsorbed on the nanotube wall. The image force acting 
on the proton moving through the nanotube embedded in a dielectric medium will 
be calculated using a two-dimensional hydrodynamic model of the nanotube atoms 
valence electrons while the surrounding dielectric medium will be described by a suitable 
dielectric function |29j . 

After outlining the basic theory used in modeling the interaction potentials between 
the proton and the nanotube and dielectric medium, we shall discuss the results of the 
proton trajectory simulations and give the concluding remarks. 

The atomic units will be used throughout the paper unless explicitly stated 
otherwise. 

2. Theory 

The system under investigation is a proton moving through an (11, 9) single- wall carbon 
nanotube embedded in a dielectric medium. The z axis coincides with the nanotube 
axis and the origin lies in its entrance transverse plane. The initial proton speed vector 
is taken to be parallel to the z axis. We assume that the nanotube is sufficiently short 
for the proton energy loss to be neglected. The initial proton speed, v, is varied between 
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3 and 8 a.u., corresponding to the energies between 0.223 and 1.59 MeV, respectively. 
The nanotube length, L, is varied between 0.1 and 0.8 fim. 

We assume that the (repulsive) interaction between the proton and nanotube atoms 
can be treated classically using the Doyle- Turner expression for the proton-nanotube 
atom interaction potential averaged axially and azimuthally [3ll [35l [36] . The resulting 
interaction potential reads 

UreAr) = E a.&?/o(fe?rc/) exp{-6,^p + {d/2)% (1) 

where Z\ = \ and = 6 are the atomic numbers of the proton and nanotube atom, 
respectively, d is the nanotube diameter, I is the nanotube atoms bond length, r is 
the distance between the proton and nanotube axis, /q designate the modified Bessel 
function of the first kind and the 0*'' order, and = (0.115,0.188,0.072,0.020) and 
bj = (0.547,0.989, 1.982,5.656) are the fitting parameters [36]. 

The dynamic polarization of the nanotube and dielectric medium by the proton is 
treated using a two-dimensional hydrodynamic model of the nanotube atoms valence 
electrons, based on a jellium-like description of the nanotube ion cores, extended to 
include the contribution of the dielectric boundary |37l [El [29]. This model includes 
the procedures of axial and azimuthal averaging, as the repulsive interaction model. 
It gives the (attractive) interaction potential between the proton and its image. Let 
us designate the proton position at time t hj r = r(t) and the electric potential at 
point R originating from the (screened) proton, perturbing the nanotube and dielectric 
boundary, by $e^t(i?, r, t). This potential is called the external potential. The image 
interaction potential at the proton position reads 

Uatr = -Y'^ind{r,t), (2) 

where ^^^^(r, t) is the electric potential at the proton position originating from the 
polarization charges induced on the nanotube and dielectric boundary by the (screened) 
proton. This potential is called the induced potential. 

Since the repulsive and attractive interaction potentials are axially symmetric, the 
initial proton speed vector is parallel to the z axis and the proton is channeled, the 
cylindrical coordinates of the proton position are r = r{t), if = and z = vt, where 
(fo is the azimuthal coordinate of the initial proton position. Let us designate the 
nanotube and dielectric boundary radii by a = d/2 and b, respectively. Then, the 
Fourier transform of the external potential at the nanotube or dielectric boundary, i.e., 
for R = a or R = b, is 

ItiZ 

^ext{R, r) = -g{R, r, m, k)5{u - kv) exp{-imipo) (3) 

where Ent = 1 is the dielectric constant of the nanotube background, m, k and uj are 
the angular oscillation mode, longitudinal wave number and angular frequency of an 
elementary excitation of the nanotube atoms valence electrons treated as an electron 
gas, and g{R,r,m, k) is the radial Green's function. 
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The Fourier transform of the induced potential at the proton position is 



^ — (XTl 

^ind{r) = -[g{a,r,m,k) + bg{b,r,m,k)^g'{b,a,m,k)] 



+ bg{b,r,m,k)m',^t{b,r), (4) 



where 



~ ^ ^extja, r) + bg{b, a, m, k)?ft^'^^^{b, r) 

+ {a/snt) [g{a, a, m, k) + bg{b, a, m, k)Ug'{b, a, m, k)] ' 

no{k'^ + w? / a^) 

X = ^ 2 ' (6) 

a{k'^ -\- / a^) -\- (5{k'^ -\- m? / 0?) — a;(a; + i7) 



and 



AT^[ent + (£nt - e^)kbI„,{kb)K'^,{kb)] ' ^ ^ 

ha is the Fourier transform of the surface polarization charge density induced at the 
nanotube, % is the response function of the polarization charge induced at the nanotube, 
no = 0.428 a.u. is the equilibrium density of the electron gas, a = ttuq, (3 = 1/4, 
and 7 — > 0+, 3? is the response function of the polarization charge induced at the 
dielectric boundary, e,^ is the dielectric function of the dielectric medium, and and 
Km designate the modified Bessel functions of the first and second kinds and the m*'* 
order, respectively; the derivatives of functions g{R,r,m,k) and ^extiR-if) a-re taken 
with respect to their first arguments. 

The induced potential at the proton position is 



ind 



'f,t)= J2_J J ^exp{i[m(po-{(jJ-kv)t]}^ind{r). (8) 



The total interaction potential between the proton and the nanotube and dielectric 
medium reads 

U (f, t) = Urep{r) - Y*md(f, t). (9) 

The spatial distributions of channeled protons in the exit transverse plane are 
generated using a computer simulation method, which includes the numerical solving of 
the proton equations of motion in the transverse plane. The rectangular coordinates of 
the initial proton position, xq and yo, are chosen randomly from a two-dimensional 
uniform distribution with condition tq — {xQ + yg)^ < a — age, where a^c = 
[97r^/(128 ^'2)]^ Oo is the nanotube atom screening radius and oq the Bohr radius. We 
take for the nanotube atoms bond length / = 0.144 nm [38] and obtain for the nanotube 
radius a = 0.689 nm. The initial number of protons is 3 141 929. 

It has been demonstrated that proton channeling in nanotubes can be analyzed 
successfully via the corresponding mapping of the impact parameter plane to the 
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scattering angle plane [SI El [IDl [HI [211 [2S] • Analogously, we analyze here the mapping 
of the entrance transverse plane to the exit transverse plane. However, since in the 
case we investigate the total interaction potential is axially symmetric, the analysis of 
this mapping can be reduced to the analysis of the mapping of the proton radial axis 
in the entrance transverse plane to the proton radial axis in the exit transverse plane. 
Further, we can take that = and analyze only the mapping of the Xq axis in the 
entrance transverse plane to the Xq axis in the exit transverse plane. The extrema in 
this mapping are the rainbow maxima or minima, and the corresponding singularities 
in the spatial distribution of channeled protons are the rainbow singularities. 

3. Results and discussion 

We shall analyze first the spatial distributions of protons channeled in the (11, 9) carbon 
nanotubes placed in vacuum. The proton speed will be f = 3 a.u. and the nanotube 
lengths L = 0.1, 0.2 and 0.3 fim. For this proton speed the spatial distributions of 
channeled protons for the nanotubes embedded in a dielectric medium do not differ 
from the corresponding spatial distributions for the nanotubes placed in vacuum. This 
part of the study represents a continuation of our previous study of the influence of the 
dynamic polarization effect on the angular distributions of protons channeled in the (11, 
9) nanotubes placed in vacuum [2^ . 

After that, we shall explore the spatial distributions of protons channeled in the 
(11, 9) nanotubes embedded in Si02. The proton speed will be f = 5 a.u. and the 
nanotube lengths L = 0.3 and 0.5 /im. For this proton speed the spatial distributions 
of channeled protons for the nanotubes embedded in a dielectric medium differ from 
the corresponding spatial distributions for the nanotubes placed in vacuum. This part 
of the study represents a continuation of our previous study of the influence of the 
image interaction potential on the angular distributions of protons channeled in the 
(11, 9) nanotubes placed in Si02 [30]. Then, we shall consider the spatial distributions 
of protons channeled in the (11, 9) nanotubes embedded in Si02, AI2O3 and Ni. The 
proton speed will be f = 8 a.u. and the nanotube length L = 0.8 a.u. 

The nanotube radius is a = 13.01 a.u. The separation between the nanotube and 
the dielectric boundaries is 3.21 a.u. Thus, the dielectric boundary radius is 6 = a + 
3.21 a.u. We describe the surrounding Si02 by the dielectric constant of 3.9 [39]. The 
dielectric responses of the surrounding AI2O3 and Ni are modeled using the methods 
described by Arista and Fuentes [21] and Kwei et al. [lO], respectively. We have found 
that in the cases of Si02 and AI2O3 a proton moving at a speed below about 3 a.u. does 
not polarize the dielectric media. This means that the dielectric medium is completely 
screened by the nanotube and does not contribute to the image force. For a proton 
speed above about 3 a.u. the polarization of the dielectric media occurs. Consequently, 
the screening of the dielectric medium by the nanotube is incomplete and it influences 
the image force. 

Figure 1(a) shows two spatial distributions of channeled protons along the x axis 
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Figure 1. (a) Two spatial distributions of channeled protons along the x axis in the 
exit transverse plane for the proton speed of 3 a.u. and the nanotube length of 0.1 
/im. The nanotube is placed in vacuum. The dashed line corresponds to the case 
without the dynamic polarization effect and the solid line to the case with it. The size 
of a bin along the x axis is 0.3 a.u. (b) Two mappings of the xq axis in the entrance 
transverse plane to the x axis in the exit transverse plane corresponding to the two 
spatial distributions. 



in the exit transverse plane for the proton speed f = 3 a.u. and the nanotube length L 
= 0.1 /im. The nanotube is placed in vacuum. The spatial distributions correspond to 
the cases without the dynamic polarization effect and with it. In each case the spatial 
distribution contains a strong central maximum and a pair of prominent peripheral 
maxima, designated by 1 in the former case and by Ij in the latter case. In the case 
with the dynamic polarization effect the central maximum is narrower and weaker than 
in the case without it. In the former case the peripheral maxima are located at x = ±8.4 
a.u., and in the latter case at a; = ±10.2 a.u. 

Figure [H^b) shows two mappings of the Xq axis in the entrance transverse plane to 
the X axis in the exit transverse plane corresponding to the two spatial distributions of 
channeled protons shown in Fig. [T](a). One can see that in each case the mapping has 
two extrema, a maximum and a minimum. The extrema in the mapping without the 
image interaction potential are designated by 1 and the extrema in the mapping with 
it by Ij. In each case the ordinates of the extrema coincide with the abscissae of the 
maxima in the corresponding spatial distribution. Since the extrema in the mappings 
are the rainbow extrema, the maxima in the spatial distributions are in fact the rainbow 
singularities. 

If the spatial distributions of channeled protons shown in Fig. [11(a) are compared 
with the corresponding angular distributions of channeled protons [24j, obtained for 
the proton speed f = 3 a.u. and the nanotube length L = 0.1 /im, it is seen that in 
the former case the rainbow effect is much more pronounced. In this case the rainbow 
maxima appear even when the image force is not included, unlike in the latter case. 
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Figure 2. (a) Two spatial distributions of channeled protons along the x axis in the 
exit transverse plane for the proton speed of 3 a.u. and the nanotube length of 0.2 
/Ltm. The nanotube is placed in vacuum. The dashed line corresponds to the case 
without the dynamic polarization effect and the solid line to the case with it. The size 
of a bin along the x axis is 0.3 a.u. (b) Two mappings of the xq axis in the entrance 
transverse plane to the x axis in the exit transverse plane corresponding to the two 
spatial distributions. 



Figure [2](a) shows two spatial distributions of channeled protons along the x axis in 
the exit transverse plane for the proton speed v = 3 a.u. and the nanotube length L = 
0.2 nm. The nanotube is placed in vacuum. The spatial distributions correspond to the 
cases without the image force and with it. In the former case the spatial distribution 
contains a central maximum and two pairs of peripheral maxima, designated by 1 and 
2, and in the latter case a central maximum and three pairs of peripheral maxima, 
designated by Ij, 2j and 3j. In the case with the image force the central maximum is 
narrower and about two times weaker while the peripheral maxima are more prominent 
than in the case without it. In the former case the peripheral maxima are located at 
X = ±7.8 and ±11.7 a.u., and in the latter case at x = ±9.6, ±10.8 and ±12.0 a.u. 

Figure EJ^b) shows two mappings of the Xq axis in the entrance transverse plane to 
the X axis in the exit transverse plane corresponding to the two spatial distributions of 
channeled protons shown in Fig. [2]^a). It is clear that in the case without the dynamic 
polarization effect the mapping has four extrema, two maximum and two minima. They 
are designated by 1 and 2. In the case with the dynamic polarization effect the mapping 
has six extrema, three maxima and three minima. They are designated by Ij, 2j and 3j. 
In each case the ordinates of the extrema coincide with the abscissae of the maxima in 
the corresponding spatial distribution. 

There are two new features in Fig. [2] when it is compared to Fig. [TJ First, each 
of the spatial distributions of channeled protons given in Fig. [2t^a) contains a pair of 
very weak maxima, designated by 2 in the case without the image force and by 3j in the 
case with it. They correspond to the very sharp maxima and minima in the mappings 
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given in Fig. [2](b), designated by 2 and 3,, which he near the nanotube wah. It is 
noteworthy that the image force does not change much the positions and shapes of the 
very weak maxima in the spatial distribution and the very sharp maxima and minima 
in the mapping. This can be attributed to the dominance of the repulsive interaction 
potential in the region near the nanotube wall. Second, the spatial distribution given 
in Fig. without the image force included contains one pair of prominent peripheral 
maxima, designated by 1, like the corresponding spatial distribution given in Fig. [T](a). 
On the other hand, the spatial distribution given in Fig. 2(a) with the image force 
included contains two pairs of prominent peripheral maxima, designated by 1, and 2j, 
unlike the spatial distribution given in Fig. [11(a) with the image force included and 
the spatial distribution given in Fig. [2]^a) without the image force included. They 
correspond to the shallow maxima and minima in the mapping given in Fig. [2](b) with 
the image force included, designated by Ij and 2j, which lie away from the nanotube 
wall. These maxima demonstrate the presence of the image force. 

Figure [3](a) shows two spatial distributions of channeled protons along the x axis 
in the exit transverse plane for the proton velocity f = 3 a.u. and the nanotube length 
L = 0.3 fim. The nanotube is placed in vacuum. The spatial distributions correspond 
to the cases without the image interaction potential and with it. In the former case the 
spatial distribution contains a central maximum and three pairs of peripheral maxima, 
designated by 1, 2 and 3, and in the latter case a central maximum and five pairs of 
peripheral maxima, designated by Ij, 2j, 3j, 4j and 5j. The maxima designated by 
2, 3, 4i and 5j are very weak. In the case with the image interaction potential the 
central maximum is much narrower and about three times weaker while the peripheral 
maxima are much more prominent than in the case without it. In the former case the 
peripheral maxima are located at x = ±7.5, ±11.1 and ±12.1 a.u., and in the latter 
case at a; = ±8.1, ±9.9, ±10.7, ±11.4 and ±12.2 a.u. 

Figure EJ^b) shows two mappings of the xq axis in the entrance transverse plane to 
the X axis in the exit transverse plane corresponding to the two spatial distributions 
of channeled protons shown in Fig. [3](a). It is evident that in the case without the 
image force the mapping has six extrema, three maxima and three minima. They are 
designated by 1, 2 and 3. In the case with the image force the mapping has 10 extrema, 
five maxima and five minima. They are designated by Ij, 2j, 3j, 4j and 5^. The maxima 
and minima designated by 2, 3, 4j and 5, are very sharp and lie near the nanotube wall. 

Figure H] shows the spatial distributions of channeled protons in the exit transverse 
plane for the proton speed v = 3 a.u. and the nanotube length L = 0.3 fim without the 
dynamic polarization effect and with it. The nanotube is placed in vacuum. The weak 
azimuthal asymmetry of the spatial distributions would be less pronounced if the initial 
number of protons were larger. 

It is obvious from Figs. [TH that for a longer nanotube the image force makes 
the spatial distributions of channeled protons richer. Figure H] shows clearly that the 
image force causes a significant increase of the proton flux near the nanotube wall, in 
a similar way as it was observed by Zhou et al. [15j. It is important to note that this 
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Figure 3. (a) Two spatial distributions of channeled protons along the x axis in the 
exit transverse plane for the proton speed of 3 a.u. and the nanotube length of 0.3 
/Ltm. The nanotube is placed in vacuum. The dashed line corresponds to the case 
without the dynamic polarization effect and the solid line to the case with it. The size 
of a bin along the x axis is 0.3 a.u. (b) Two mappings of the Xq axis in the entrance 
transverse plane to the x axis in the exit transverse plane corresponding to the two 
spatial distributions. 



increase of the proton flux occurs at tlie distances from the nanotube wall of the order 
of 0.2 nm, coinciding with the typical separations between the nanotube wall and the 
atoms that can be adsorbed on them. Therefore, it is evident that the image force 
makes possible the use of proton channeling for detecting and locating the atoms and 
molecules adsorbed on the nanotube wall. 

Figure [5] gives five proton trajectories in the xz plane for the Xq coordinates of the 
initial proton position corresponding to the five maxima in the spatial distribution of 
channeled protons given in Fig. [3](a) with the image interaction potential. The analysis 
shows that these trajectories, which are called the rainbow trajectories, can be classified 
according to the number of proton deflections within the total interaction potential well 
|24j . The number of proton deflections determines the order of the rainbow. One can 
see that trajectory 4j, corresponding to maximum 4j in the spatial distribution, includes 
one proton deflection. Therefore, the rainbow it represents is the primary rainbow. 
The rainbows represented by trajectories !«, 3j and 5i, corresponding to maxima Ij, 
3j and 5j in the spatial distribution, are the secondary rainbows, and that represented 
by trajectory 2j, corresponding to maxima 2j in the spatial distribution, is the tertiary 
rainbow. 

We next compare the effects of the image force on the spatial distributions of protons 
channeled in the nanotube placed in vacuum and embedded in Si02. The proton speed 
is increased to f = 5 a.u. and the nanotube lengths are L = 0.3 and 0.5 fim. 

Figure [6](a) shows two spatial distributions of channeled protons along the x axis in 
the exit transverse plane for the proton speed v = 5 a.u. and the nanotube length L = 
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Figure 4. Two spatial distributions of channeled protons in the exit transverse plane 
for the proton speed of 3 a.u. and the nanotube length of 0.3 /im. The nanotube is 
placed in vacuum. The spatial distributions correspond to the cases (a) without the 
dynamic polarization effect and (b) with it. The sizes of a bin along the x and y axes 
are 0.3 a.u. 
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Figure 5. Five proton trajectories in the xz plane for the xq coordinates of the 
initial proton positions corresponding to the five maxima in the spatial distribution of 
channeled protons shown in Fig. 3(a) with the dynamic polarization effect included. 

0.3 /iin. They correspond to the cases in which the nanotube is placed in vacuum and 
embedded in Si02. The image force is included. In each case the spatial distribution 
contains a central maximum and two pairs of peripheral maxima, designated by Ij and 
2j in the former case, and by l^^ and 2^^ in the latter case. The maxima designated by 
2j and 2^ are very weak. In the case in which the nanotube is placed in vacuum the 
peripheral maxima are located at x = ±9.3 and ±11.9 a.u., and in the case in which it 
is embedded in Si02 at x = ±8.7 and ±11.9 a.u. 

Figure Mjo) shows two mappings of the Xq axis in the entrance transverse plane to 
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Figure 6. (a) Two spatial distributions of channeled protons along the x axis in the 
exit transverse plane for the proton speed of 5 a.u. and the nanotube length of 0.3 fiui. 
The dashed line corresponds to the case in which the nanotube is placed in vacuum 
and the solid line to the case in which it is embedded in Si02- The size of a bin along 
the X axis is 0.3 a.u. (b) Two mappings of the Xq axis in the entrance transverse plane 
to the X axis in the exit transverse plane corresponding to the two spatial distributions. 



the X axis in the exit transverse plane corresponding to the two spatial distributions of 
channeled protons shown in Fig. [6](a). One can see that in each case the mapping has 
four extrema, two maxima and two minima. The extrema in the mapping in the case 
in which the nanotube is placed in vacuum are designated by Ij and 2j, and in the case 
in which it is embedded in Si02 by 1^ and 2^. The maxima and minima designated 
by 2j and 2^ are very sharp and lie near the nanotube wall. In each case the ordinates 
of the extrema coincide with the abscissae of the maxima in the corresponding spatial 
distribution. Thus, the maxima in the spatial distributions are in fact the rainbow 
singularities. 

Figure [7|(a) shows two spatial distributions of channeled protons along the x axis in 
the exit transverse plane for the proton speed v =5 a.u. and the nanotube length L = 
0.5 yum. They correspond to the cases in which the nanotube is placed in vacuum and 
embedded in Si02. The dynamic polarization effect is included. In the former case the 
spatial distribution contains a central maximum and four pairs of peripheral maxima, 
designated by Ij, 2j, 3j and 4j, and in the latter case a central maximum and three pairs 
of peripheral maxima, designated by 1^, 2^ and 3^. The maxima designated by 3j, 4j, 
2^ and 3^ are very weak. In the case in which the nanotube is placed in vacuum the 
peripheral maxima are located at x = ±9.0, ±9.9, ±11.2 and ±12.1 a.u., and in the 
case in which it is embedded in Si02 at x = ±9.3, ±11.2 and ±12.1 a.u. 

Figure ID^b) shows two mappings of the Xq axis in the entrance transverse plane to 
the X axis in the exit transverse plane corresponding to the two spatial distributions 
of channeled protons shown in Fig. [7](a). It is evident that in the case in which the 
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Figure 7. (a) Two spatial distributions of channeled protons along the x axis in the 
exit transverse plane for the proton speed of 5 a.u. and the nanotube length of 0.5 fim. 
The dashed line corresponds to the case in which the nanotube is placed in vacuum 
and the solid line to the case in which it is embedded in Si02. The size of a bin along 
the X axis is 0.3 a.u. (b) Two mappings of the xq axis in the entrance transverse plane 
to the X axis in the exit transverse plane corresponding to the two spatial distributions 
of channeled protons. 



nanotube is placed in vacuum the mapping has eight extrema, four maxima and four 
minima. They are designated by Ij, 2j, 3, and 4j. In the case in which the nanotube 
is embedded in Si02 the mapping has six extrema, three maxima and three minima. 
They are designated by 1^, 2^ and 3^. The maxima and minima designated by 3j, 4^, 
2rf and 3^ are very sharp and he near the nanotube wall. In each case the ordinates 
of the extrema coincide with the abscissae of the maxima in the corresponding spatial 
distribution. 

It is interesting to note that the spatial distribution of channeled protons given in 
Fig. Wl^a.) when nanotube is embedded in Si02 contains one pair of prominent peripheral 
maxima rather than two such pairs of maxima, as in the other spatial distribution given 
in this figure (when nanotube is placed in vacuum). This is explained by the weakening 
of the image force at this proton speed due to the presence of the dielectric medium, as 
it was shown by Borka et al. [30] . 

Let us now compare the spatial distributions of channeled protons along the x axis 
in the exit transverse plane for the proton speed v = 3 a.u. and the nanotube length L 
= 0.3 fim, for v = 5 a.u. and L = 0.5 fim, and for f = 8 a.u. and L = 0.8 fim. The 
dynamic polarization effect is included. These spatial distributions are characterized by 
the same duration of the process of proton channeling, i.e., by the same proton dwell 
time. In the case in which v = 3 a.u. and L = 0.3 /xm the spatial distribution with 
the nanotube placed in vacuum contains five rainbow maxima [see Fig. [3](a)], as the 
spatial distribution with the nanotube placed in Si02. In the case in which v = 5 a.u. 
and L = 0.5 fim the spatial distribution with the nanotube placed in vacuum contains 
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Figure 8. (a) Four spatial distributions of channeled protons along the x axis in 
the exit transverse plane for the proton speed of 8 a.u. and the nanotube length of 
0.8 Ilia. The dashed line corresponds to the case in which the nanotube is placed in 
vacuum, and the solid, dotted and dash-dotted lines correspond to the cases in which 
the nanotube is embedded in Si02, AI2O3 and Ni, respectively. The size of a bin along 
the X axis is 0.3 a.u. (b) Four mappings of the xq axis in the entrance transverse plane 
to the X axis in the exit transverse plane corresponding to the four spatial distributions. 



four rainbow maxima, and the spatial distribution with the nanotube placed in Si02 
three rainbow maxima [see Fig. H^a)]. In the case in which v = 8 a.u. and L = 0.8 
Hm the spatial distribution with the nanotube placed in vacuum contains three rainbow 
maxima, as the spatial distribution with the nanotube placed in Si02 [see Fig. EJ^a)]. 
It is clear that for the same proton dwell time the number of rainbow maxima in the 
spatial distribution decreases with the proton speed. 

Figure [H](a) shows four spatial distributions of channeled protons along the x axis 
in the exit transverse plane for the proton speed v = 8 a.u. and the nanotube length 
L = 0.8 fim. They correspond to the cases in which the nanotube is placed in vacuum 
and embedded in Si02, AI2O3 and Ni. The image interaction potential is included. In 
each case the spatial distribution contains a strong central maximum and three pairs 
of peripheral maxima, designated by Ij, 2j and 3j in the case in which the nanotube is 
placed in vacuum, and by 1^, 2^ and 3^ in the cases in which it is embedded in Si02, 
AI2O3 and Ni. It is clear that the spatial distributions in the cases in which the nanotube 
is embedded in the three dielectric media are very similar to each other, and that they 
do not differ significantly from the spatial distribution in the case in which it is placed 
in vacuum. The spatial distributions in the cases in which the nanotube is embedded 
in Si02 and Ni almost coincide. 

Figure Mjo) shows four mappings of the xq axis in the entrance transverse plane to 
the X axis in the exit transverse plane corresponding to the four spatial distributions of 
channeled protons shown in Fig. [H](a). One can see that in each case the mapping has 
six extrema, three maxima and three minima. The extrema in the mapping in the case 
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Figure 9. Eight proton trajectories in the x{z/v) plane for the xq coordinates of the 
initial proton positions of 3 and 8 a.u. and the proton speeds of (a) 5 a.u. and (b) 
8 a.u. The dashed lines correspond to the cases in which the nanotube is placed in 
vacuum and the solid lines to the cases in which it is embedded in Si02. 



in which the nanotube is placed in vacuum are designated by 1, 2j and Sj, and in the 
cases in which it is embedded in Si02, AI2O3 and Ni by 1^, 2^ and 3^. 

Let us now consider two different types of proton trajectories in the channehng 
process under consideration. One should be aware of the facts that in the case without 
the attractive interaction potential the protons can oscillate transversely only between 
the sides of the nanotube wall, while in the case with it they can oscillate transversely 
between the sides of the nanotube wall and also around the two minima of the total 
interaction potential, which appear due to the presence of the attractive interaction 
potential. 

Figure [9] shows the proton trajectories in the x{z/v) plane for the xq coordinates 
of the initial proton positions of 3 and 8 a.u. and the proton speeds f = 5 and 8 a.u. 
They correspond to the cases in which the nanotube is placed in vacuum and embedded 
in Si02. The image force is included. The values of |xo| are sufficiently small, i.e., 
smaller than about 11 a.u., to prevent the transverse proton oscillations between the 
sides of the nanotube wall. Instead, the protons oscillate transversely around the two 
minima appearing due to the presence of the attractive interaction potential, having 
the X coordinates about ±9.3 a.u. It is obvious that the amplitudes of the transverse 
proton oscillations are smaller and their periods larger for the larger proton speed {v 
= 8 a.u.). This is attributed to the weakening of the attractive interaction potential 
with the proton speed. It is also clear that for these values of |a;o| the influence of the 
dielectric medium on the proton trajectories is significant. For the smaller proton speed 
(f = 5 a.u.) the presence of the dielectric medium causes an increase of the period of 
transverse proton oscillations while for the larger proton speed (f = 8 a.u.) the situation 
is opposite. 
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Figure 10. Four proton trajectories in the x{z/v) plane for the xq coordinate of the 
initial proton position of 12 a.u. and the proton speeds of (a) 5 a.u. and (b) 8 a.u. 
The dashed lines correspond to the cases in which the nanotube is placed in vacuum 
and the solid lines to the cases in which it is embedded in Si02. 



Figures [TOT a) and (b) show the proton trajectories in the x{z/v) plane for the 
Xq coordinate of the initial proton position of 12 a.u. and the proton speeds v = 5 
and 8 a.u. They correspond to the cases in which the nanotube is placed in vacuum 
and embedded in Si02. The dynamic polarization effect is included. The value of 
|xo| is sufficiently large, i.e., larger than about 11 a.u., to enable the transverse proton 
oscillations between the sides of the nanotube wall. Since the nanotube wall is defined 
by the repulsive interaction potential of the Doyle- Turner type, the proton trajectories 
resemble those in the case of a box with rigid walls. It is also evident that for this value 
of |xo| the infiuence of the dielectric medium on the proton trajectories is very small, 
even after a number of transverse proton oscillations between the nanotube walls. 

4. Conclusions 

We have presented theoretical investigations of the effects of dynamic polarization of 
the nanotube atoms valence electrons on the spatial distributions of protons channeled 
through the (11, 9) single-wall carbon nanotubes placed in vacuum and embedded in 
various dielectric media, for the proton speeds between 3 and 8 a.u. and the nanotube 
lengths between 0.1 and 0.8 fim. This study complements our previous investigations 
of the effects of dynamic polarization on the angular distributions of protons channeled 
through the nanotubes, where we found that the image force gives rise to the rainbow 
effect [Ml [25] , which can be strongly affected by the surrounding dielectric medium for 
the proton speeds above 3 a.u. [30j. Similarly, our present study has revealed that the 
prominent rainbow maxima exist in the spatial distributions of channeled protons, which 
correspond to the well-defined extrema in the mapping of the entrance transverse plane 
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to the exit transverse plane. While the rainbows in the angular distributions of channeled 
ions can be easily measured and applied to provide information on the structure and 
atomic forces inside the nanotubes, the rainbows in the spatial distributions of channeled 
ions can be employed for detecting and locating the atoms and molecules intercalated 
in the nanotubes. 

Specifically, besides the strong central maxima in all the spatial distributions of 
channeled protons, we have found that, as the nanotube length increases, the very weak 
maxima lying near the nanotube wall are present too, without and with the image force 
and the dielectric media included. However, our most important finding is the prominent 
peripheral maxima in the spatial distributions at the distances from the nanotube wall 
of the order of a few tenths of a nanometer. These distances are perfectly suitable for 
applying proton channeling to probe the atoms and molecules adsorbed on the nanotube 
wall. We have revealed that the image force is responsible for the appearance of the 
additional prominent peripheral maxima as well as for making them more prominent at 
the expense of the central maximum, as the nanotube length increases. 

Since the spatial distribution of channeled protons gives us a detailed information 
about the proton flux within the nanotube, it appears that a careful studying of the 
speed dependence of the image force as well as of the properties of the nanotube and 
its surrounding can help us better understand and perhaps even find a way to induce 
a spatial redistribution of channeled protons towards the nanotube wall to be used for 
probing the atoms and molecules intercalated in the nanotube. 

Besides such diagnostic apphcations, we would like to mention the possibility of 
producing nanosized ion beams with the nanotubes embedded in various dielectric media 
for applications in biomedical research [6l [12], [13]. Our findings can shed more light 
on the issue of shape of the spatial distribution of channeled ions suitable for such 
applications, especially if there are concerns about the increase of the proton flux in the 
peripheral region of a nanosized ion beam at the expense of the flux in its central region. 
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